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ABSTRACT 

We present three Suzaku observations of the black hole candidate 4U 1957+1 1 (V1408 Aql) — a source 
that exhibits some of the simplest and cleanest examples of soft, disk-dominated spectra. 4U 1957+1 1 also 
presents among the highest peak temperatures found from disk-dominated spectra. Such temperatures may be 
associated with rapid black hole spin. The 4U 1957+1 1 spectra also require a very low normalization, which 
can be explained by a combination of small inner disk radius and a large distance (> lOkpc) which places 
4U 1957+11 well into the Galactic halo. We perform joint fits to the Suzaku spectra with both relativistic 
and Comptonized disk models. Assuming a low mass black hole and the nearest distance (3Mq, lOkpc), the 
dimensionless spin parameter a* = Jc/ GM^ > 0.9. Higher masses and farther distances yield a* Ri 1. Similar 
conclusions are reached with Comptonization models; they imply a combination of small inner disk radii (or, 
equivalently, rapid spin) and large distance. Low spin cannot be recovered unless 4U 1957+1 1 is a low mass 
black hole that is at the unusually large distance of > 40kpc. We speculate whether the suggested maximal 
spin is related to how the system came to reside in the halo. 

Subject headings: accretion, accretion disks - black hole physics - radiation mechanisms:thermal - X- 
rays:binaries 


1. INTRODUCTION 

4U 1957+11 is one of the few black hole candidates 
(BHC), like LMC X-L LMC X-3, and Cyg X-1 , that histor- 
ically has been persistently active. In contrast to these other 
sources, however, it is a Low Mass X-ray Binary (LMXB; 
see Thorstensen 1987; Hakala et al. 1999; Russell et al. 2010; 
Bayless et al. 2011). Also in contrast to other persistent 
Galactic BHC sources, 4U 1957+1 1 apparently has remained 
always in a spectrally soft state (Nowak et al. 2008). Its soft 
spectrum is well-modeled as a simple disk (Nowak & Wilms 
1999; Wijnands et al. 2002; Nowak et al. 2008; Dunn et al. 
2010), i.e., a multi-temperature blackbody (Mitsuda et al. 
1984) characterized by a peak temperature and a normaliza- 
tion related to the disk inner radius, inclination, and object 
distance. Approximately 15% of the over 40 pointed obser- 
vations taken with the Rossi X-ray Timing Explorer (RXTE), 
however, also show a significant contribution (20-75% of the 
observed 3-18 keV flux) from a non-thermal/steep powerlaw 
component (Nowak et al. 2008; Dunn et al. 2010). 

Optical observations suggest that we may be viewing the 
disk in 4U 1957+1 1 at a high inclination of ~ 75°. Modula- 
tion over a 9.33 hr orbital period has ranged from =c 1 0% and 
sinusoidal (Thorstensen 1 987) to ±30% and complex (Hakala 
et al. 1999). The latter has been interpreted as a high inclina- 
tion warped disk being partly occulted by the secondary. On 
the other hand, Russell et al. (2010) failed to find any opti- 
cal modulation over the orbital period, while Bayless et al. 
(201 1) found nearly sinusoidal modulation. These latter au- 
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thors attributed the optical modulation to illumination of 
the secondary, and noting a lack of any flattening at the 
lightcurve minimum, they modeled the system inclination 
with values as low as 20° and not greater than 70° . 

Currently, the mass, distance, and inclination of 
4U 1957+11 are unknown. Based upon a comparison 
of the optical flux to the estimated optical luminosities of 
other BHC, however, Russell et al. (2010) have argued that 
4U 1957+1 1 must lie at a distance > 20kpc. The possibility 
of a high inclination, however, makes such comparisons 
problematic. Russell et al. (2011) recently have presented 
near-simultaneous radio {Expanded Very Large Array - 
EVLA), optical (Eaulkes Telescope North), and UV/X-ray 
(Swift) observations of 4U 1957+11. These observations 
show the soft state flux to be extremely quenched (factors 
of 330-810 times lower than the expected radio/X-ray flux 
ratio for hard state black hole sources), and demonstrate 
that the optical/UV flux are roughly consistent with the 
expectation of an irradiated disk. Those observations, 
however, do not definitively constrain the source’s inclina- 
tion or distance. A minimum distance of > 5 kpc; however, 
has been determined via high resolution X-ray spectroscopy. 
Neix 13.45 A absorption, associated with the warm/hot 
phase of the interstellar medium, is detected with a sufficient 
equivalent width to suggest that 4U 1957+11 resides above 
the galactic plane at this minimum distance or beyond (Yao 
et al. 2008; Nowak et al. 2008). 

Nowak et al. (2008) showed that, with the exception of the 
w 15% of the time when a non-thermal component becomes 
prominent, disk models of the soft X-ray component have 
an approximately constant radius but a temperature that in- 
creases with flux (see Fig. 9 of Nowak et al. 2008). There 
is some indication for a slight increase in disk radius at the 
lowest flux levels, which was interpreted as possibly being 
associated w'ith the source approaching a transition to a spec- 
trally hard state. Such a transition is expected to occur at 
2-3%LEdd, where LEdd is the Eddington luminosity for the 
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source in question (Maccarone 2003). At the other extremes 
of flux in 4U 1957+1 1 , the simple disk spectrum “drops out” 
(disk radius decreases in disk+powerlaw fits, or disk temper- 
ature decreases in disk+Compton corona fits) and is replaced 
by the non-thermal/steep powerlaw component (Nowak et al. 
2008). 

Although other persistent systems transiently show disk- 
dominated spectra' , the persistent, soft spectra of 4U 1957+1 1 
are perhaps the simplest and cleanest examples of “disk spec- 
tra”. XMM-Newton and Chandra observations show that 
there is very little absorption (Ah = 1-2 x 10^' cm"^), while 
RXTE observations show that a hard tail contributes usually 
< 20% of the flux. Furthermore, any additional components 
such as a broad Fe line (typically < 2% peak residud in the 
fits) or a “smeared Fe edge” (which was never required or in- 
dicated in any of the RXTE spectra; Nowak et al. 2008) are 
very weak or absent. Thus the soft spectrum of 4U 1957+1 1 
becomes the ideal testbed for modern disk atmosphere mod- 
els that incorporate spin and other General Relativistic effects 
into their calculations (e.g., Li et al. 2005; Davis et al. 2005, 
2006; Shafee et al. 2006). 

Nowak et al. (2008) showed that the disk-fits to the XMM- 
Newton, Chandra, and RXTE spectra of 4U 1957+11 are 
characterized by a very low normalization, indicating some 
combination of large distance and low compact object mass, 
and very high inner disk temperature (1.3-1.8keV). In other 
sources the latter has been associated with high black hole 
spin. To assess this possibility, Nowak et al. (2008) applied 
the relativistic disk model, kerrbb, of Li et al. (2005). The 
model parameters include a compact object mass and dis- 
tance, disk inclination, accretion rate, spectral hardening fac- 
tor (the ratio of color temperature to effective temperature), 
and dimensionless spin of the black hole, a* = Jc/GM^. For 
a wide variety of masses and distances, fits to the XMM- 
Newton and Chandra spectra preferred maximal spin, a* « 1 . 
The RXTE spectra, in fact, could not be described with low 
spin, even when including a coronal or powerlaw component. 

Various theoretical models of a disk atmosphere have sug- 
gested a spectral hardening factor of ~ 1.7. If this is correct, 
and utilizing the fact that the lowest flux observations must 
be ^ 3%LEdd, then the spectral fits of 4U 1957+1 1 strongly 
suggest that it is a maximally spinning, w 1 6 Mq black hole at 
a distance of «22kpc (Nowak et al. 2008). In this work we 
explore that claim with three recent Suzaku observations. 

In Section 2 we describe the Suzaku observations and the 
details of the data processing. We then present the model 
fits to the resulting spectra in Section 3 . Section 4 describes 
the scaling relations that we used to translate fits using a low 
mass, short distance black hole to comparably good fits in- 
volving a larger mass, further distance black hole. We then 
present our conclusions in Section 5 . 

2. OBSERVATIONS AND DATA ANALYSIS 

The Suzaku data were reduced with tools from the HEA- 
SOFT v6.9 package and the calibration files current as of 20 1 0 
September. The instruments on Suzaku (Mitsuda et al. 2007) 
are the X-ray Imaging Spectrometer (XIS; Koyama et al. 
2007) CCD detector covering the « 0.3-1 0 keV band, and the 
Hard X-ray Detector (HXD; Takahashi et al. 2007) comprised 
of the PIN diode detector (PIN) covering the w 10-70keV 

' LMC X-1 is likely wind-fed, and, for a soft state, shows an unusual and 
highly variable disk spectrum (Now'ak et al. 2001). LMC X-3 cycles through 
soft/hard state transitions {Wilms et al. 2001). 


TABLE 1 

LOGOF4U 1957+11 Observations 


Date 

(yyyy-mm-dd) 

(MJD) 

ObsID 

Exposure 

(ksec) 

2010-05-04 

55320.82 

405057010 

35.80 

2010-05-17 

55333.80 

405057020 

37.20 

2010-11-01 

55502.19 

405057030 

15.50 


Note. — Exposure times are after good time filter- 
ing, and represent the maximum exposure in an XIS de- 
tector when combining all data modes. MJD is the (non- 
barycenter corrected) average Modified Julian Date of 
the observation. 



Fig. 1 .— The aspect corrected XIS 1 0.5-8 kev image on a logarithmic 
intensity scale from ObsID 405057010. The boxes outline the detector re- 
gion from which the spectra were generated. The length of the short edge of 
the outer box is approximately 4'. 5. The inner boxed regions with diagonal 
slashes were excluded from the spectra. 

band and the gadolinium silicate crystal detector (GSO) cov- 
ering the w 60-600 keV band. The XIS had four separate 
detectors, XIS 0-3, with XIS 1 being a backside illuminated 
CCD. XIS 2 was lost due to a micrometeor hit in late 2006, 
and thus we only consider data from the remaining three XIS 
detectors. For all three observations described here, the hard 
X-ray fluxes were too low to lead to a significant detection in 
either of the HXD detectors, and therefore we do not include 
these data. A log of the observations is found in Table I . 

All of the XIS observations were run in a sub-array mode 
wherein 1/4 of the CCD was exposed with 2 second inte- 
gration times. In preparing the XIS spectra, we first cor- 
rected each detector for Charge Transfer Inefficiency (CTl) 
using the xispi tool, and then reprocessed the data with 
xselect using the standard xisrepro selection criteria. 
Due to thermal flexing of the spacecraft, the attitude of the 
Suzaku spacecraft exhibits variability over the course of the 
observations and therefore the image of the source is not at a 
fi.xed position on the CCD. Standard processing reduces this 
variability and improves the Point Spread Function (PSF) im- 
age (Uchiyamaet al. 2008); however, the standard tools do not 
yet fully correct the variable attitude induced blurring. We 
created an improved attitude reconstruction, and hence nar- 
rower PSF images, using the aeattcor.sl software de- 
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FiG.2.— The combined XIS 0, 1, and 3 spectra from each of the three 
Suzaku observations of 4U i 957+11, flux corrected with the detector re- 
sponses. (ObsID 405057010 is the top spectrum, obsID 405057020 is the 
middle spectrum, and obsID 405057030 is the bottom spectrum.) The three 
observations share a common binning of approximately half width half max- 
imum resolution of the detectors. Lighter colored data points represent 
energy bands that were ignored in the spectral fits described below. 

scribed^ by Nowak et al. (201 1). For the third observation, 
however, approximately half of the good time intervals (i.e., 
half of each spacecraft orbit) were affected by wobbles severe 
enough (> 1') such that a large fraction of the image on the 
XIS 1 detector fell beyond the boundaries of the detector’s 
active regions. Although the images did not fall beyond the 
active region boundaries on the XIS 0 and XIS 3 detectors, 
there are serious questions about the accuracy of the gener- 
ated response functions for such large off axis angles. We 
therefore excluded these times from the spectra, and reflect 
this in Table 1 . 

Using the image generated from the improved attitude cor- 
rection, we then used the pile_estiniate . si S-Lang 
script^ (Nowak et al. 201 1) to estimate the degree of pileup in 
the XIS observations. For the spectra described in this work, 
the center of the PSF images were affected by as much as a 
33% pileup fraction. We therefore exclude two, overlapping 
rectangular regions (« 15 x 60 pixels, nearly at rightangles to 
one another) aligned along the brightest portions of the cen- 
tral PSF (see Fig. 1). We estimate that the remaining residual 
pileup fractions are <3%. 

The outer limit of the extraction area was a rectangular re- 
gion limited by the 4' width of the 1/4 array boundaries. Ad- 
ditionally, we extracted background spectra from small rect- 
angular regions, near the edges of the chips, that avoided as 
much as possible both the 4U 1957+1 1 image as well as the 
detector images of the on-board calibration sources. It is not 
entirely clear to what degree these regions are dominated by 
“true background”, as opposed to the wings of the PSF from 
the source. However, the normalized background flux remains 
at < 2% of the source flux between 0.5-8 keV, and thus does 
not strongly affect the analyses presented in this work. The 
background flux starts to become an increasingly larger frac- 
tion of the source flux above 8 keV. For this reason, as well as 
due to any other uncertainties in the high energy calibration 
of the XIS detectors, we do not consider data above 8 keV. 

Events in the XIS detectors were read out from either 3x3 
or 5 X 5 pixel islands. We created individual spectra and re- 
sponse files for each detector and data mode combination. Lit- 

- See also: http://space.mit.edu/CXC/software/suzaku/ 
aeatt.html- 

http://space.mit.edu/CXC/software/suzaku/pest.html 
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Fig. 3 Ratios of the flux corrected spectra shown in Fig. 2. The top ratio 
is observation 1/observation 3, the middle line is observation l/observation2, 
and the bottom line is observation 2/observation 3. 

tie source variability was detected over the course of each ob- 
servation, therefore the spectra were created for the full good 
time intervals represented by standard processing (excluding 
the times of very bad attitude errors, discussed above). Re- 
sponse matrices and effective area files were created with the 
xisrmf gen and xissimarf gen tools, respectively. For 
each observation , spectra for the different readout modes were 
combined during analysis on a chip-by-chip basis using the 
combine_datasets function within the Interactive 
Spectral Interpretation System (ISIS) fitting 
program (Houck & Denicola 2000). That is, observations 
were kept separate, but within an observation all XIS 0 spec- 
tra were combined together, all XIS 1 spectra were combined, 
etc. 

Although spectra for each observation and each XIS chip 
were kept separate, all spectra were jointly grouped on a com- 
mon grid such that the XIS 0 spectra for the third (i.e., the 
faintest and shortest integration time) observation had a min- 
imum combined signal-to-noise ratio of 5 in each energy bin 
(w 25 total counts in each bin, although the estimated back- 
ground was included in the S/N criterion) and that the mini- 
mum number of channels per energy bin was at least the half 
width half maximum of the spectral resolution (see Nowak 
et al. 201 1 ,for a further explanation of this criterion). 

In Fig. 2 we present the flux corrected spectra, i.e., the 
counts spectra divided by the integrated response (see the 
ISIS user manual for a full description of this process). The 
spectra show strong deviations from a smooth continuum in 
the 1 .5-2.6 keV region, near X-ray spectral features of Si and 
Ir. These systematic effects are seen in a wide variety of un- 
related Suzaku spectra, hence our decision to exclude them 
from our fits. Likewise, there is a strong spectral feature near 
0.55 keV, which we also suspect of being systematic in na- 
ture, hence our exclusion of data below 0.7 keV. To avoid 
these regions of poorly understood response, we only consid- 
ered spectral energy ranges 0.7-1 .5 keV and 2.6-8 keV. The 
smaller features near 0.75 keV and 0.9 keV are plausibly as- 
sociated with Fe L, Ne K, and Ne IX absorption in the inter- 
stellar medium (ISM), as we discuss below. 

The suspected systematic features are relatively constant 
among the three observations, as is evident from the ratio of 
the flux corrected spectra as shown in Fig. 3. None of the 
systematic features are seen in these ratios. The ratios are rel- 
atively flat < 2 keV, and show a curve upward, with the degree 
of this curving being greatest for the ratio of the flux corrected 
spectra for the first and third observation. We discuss this fur- 
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Fig. 4.— Expectations for the flux ratio in a diskpn model when solely 
changing the peak temperature of the disk. The reference model has a 1 3 keV 
peak temperature. The peak temperature is varied in steps of 0.05 keV, up to 
1.5keV. 



Biergv‘ (ke\') 

Fig. 5.— The joint fit of an eqpair Comptonizalion model to the three 
spectra. The spectra share a common interstellar absorption column and a 
common normalization (i.e., disk inner radius), but have independent disk 
temperatures and coronal parameters. Black/grey spectra are the combined 
XfS 0 and 3 spectra, and blue/light blue are the XIS 1 spectra. Lightly shaded 
data were not included in the fit. The top residuals panel shows the ratio resid- 
uals for obsID 405057010, the middle residuals panel shows the ratio residu- 
als for obsID 405057020, and the bottom panel shows the ratio residuals for 
obsID 405057030. 

ther in the next section. 

3. SPECTRAL FITS 

The spectral shape shown in Fig. 2, and the ratios seen 
in Fig. 3, are those expected for a disk-dominated spectrum 
wherein the flux changes are primarily dominated by varia- 
tions of the peak temperature of the disk. We show the theo- 
retical expectation for such a disk in Fig. 4. Here we use the 
diskpn model (Gierlihski et al. 1999), which essentially is 
the spectrum from a Shakura-Sunyaev type disk (Shakura & 
Sunyaev 1973) calculated using a simple pseudo-Netwonian 
potential and a no-torque boundary condition at the disk inner 
edge located at (yGM/c^. The flux goes to zero at the disk in- 
ner edge, and the disk temperature peaks at « lOGMc^. This 
basic disk spectrum will underly the Comptonizalion fits dis- 
cussed below. 

At energies approximately below the peak temperature. 



Energy (keV) 

Fig. 6.— A simultaneous fit to all three observations using the eqpair 
model, but excluding the broad Gaussian line in the Fe region. (The spectra 
were combined for the plot, but not the fits. 

dominated by the Rayleigh-Jeans part of the thermal spec- 
tra, the ratios are fairly flat. At energies above the disk peak 
temperature, on the Wien tail of the thermal spectra, the ratios 
curve upward. There is a very good correspondence between 
the curves in Figs. 3 and 4. However, based upon our ex- 
perience with RXTE spectra which extended out to w 20 keV, 
there is the possibility of additional spectral hardening beyond 
that associated with any rise in disk temperature. We therefore 
consider the eqpair Comptonizalion model (Coppi 1999) 
which can use the diskpn as its input seed photon spectrum. 

Given that the three spectra appear commensurate with a 
constant disk radius, and that such a nearly constant disk 
radius is consistent with our experience for « 85% of the 
RXTE spectra that we have previously examined (see Fig. 9 
of Nowak et al. 2008), we explore a model wherein we fit 
all three spectra simultaneously with an eqpair model for 
which we constrain the eqpair normalization (i.e., the inner 
disk radius) to be the same for each spectrum. The physical 
implications of this common normalization are discussed 
in §5. 

In addition to constraining the eqpair normalization, 
we also constrain the absorption to be the same for each ob- 
servation, which we here model with the TBnew model (an 
updated version of the model of Wilms et al. 2000, using the 
interstellar abundances from that work) and using a fixed Ne 
IX 13 .447 A absorption line with equivalent width of 7 A. The 
latter is consistent with the high spectral resolution studies of 
Yao et al. (2008) and Nowak et al. (2008) (see the below dis- 
cussion of the ISM absorption). Excluding this component 
increases the fitted by « 12. 

For each observation we allow the disk peak temperature, 
^Tdisk, the coronal compactness, 4/4. and the coronal seed 
optical depth^, r, , to be independent parameters. We include 
Gaussian lines to the fits; however, as their presence is not 
strongly confirmed and there is no good evidence for a nar- 
row line in either these spectra nor in the previously exam- 
ined Chandra-HETG spectra (Nowak et al. 2008), we tie all 
the line energies (constrained to lie between 6.2-6.9keV) and 
widths (fixed to 0.3 keV), and let only the line normalizations 

Table 2, however, presents the nei optical depth, as eqpair includes 
calculation of the equilibrium pair production in the corona. 
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TABLE 2 

constant*TBnew* (eqpair+gaussian) Fits TO 4U 1957+11 Spectra 


ObslD 

Nu 

Mrqp 

^'Tdisk 

4/4 

7"nel 

At-y 

Qiso 

0kis3 

0.7-8 keV Hux 

xVdof 


(10^* cm -) 

(10-*) 

(keV) 



(10-*) 



(lO'^ergcm-^s) 


405057010 

405057020 

405057030 

U.-.UU_^QQJ 

1 Q^5-H).00i 

l.y-0_Qooi 

1 477+<).00I 
'-0.00I 

1 31A+O.OOI 

Q 37+0.03 
'-0.02 

0 05^-^^ 
^•^^-0.01 

o.o4!S:gl 

'-O.002 
•7 4+0-5 
“•^-0.1 
4.4*:| 

0.1/04 

1 01 

‘■^^^-0.002 
1 n7S+0-l^2 

1 a^7+0.002 
^ '-0.002 

1 

i •^-'4 _(),o() 2 
1 ns^+0.002 

i .KJD /_<) Q05 

1.09 + 0.01 
0.80 ±0.04 
0.72 ±0.03 

1750.4/1014 

405057010 

405057020 

405057030 

n iqq+O-007 
’-'•^^-0.008 

1 926'^-^' 

t 47z:+0.0()l 

1 343**^-^^ 
1 310'^^-^^ 

0.38^:“ 

0 05^ *^^ 
^•^^-0,01 
A Ar+0,01 
^•^^-0.01 

o.oirji:^ 

2<:l 

4.5^:1 

1 4+0.4 
-Q-4 

^ •^‘*'-'_0.002 
1 077^-^“ 

1 fV^'+0-003 
i.UOO_^(y^ 

i ,uzu_Q 
i.Uoi_^ooi 

1 A«^+0.004 


1016.9/1014 


Note. — Errors are 90% confidence level for one interesting parameter, except for the 0.8-S keV (absorbed) flux error which is from the standard deviation 
of the XIS detector cross normalizations. Nn is the neutral column; Afcqp is the eqpair normalization (see text); A'Tdjsk is the peak temperature of the disk seed 
photons; 4/4 is the relative corona to disk compactness {is is fixed to 1); Tum is the net coronal optical depth (including pair production); N-y is the Gaussian 
line normalization in units of photons cm“- s“* ; and CjIsk are the normalization constants for the XIS 0 and 3 detectors, relative to the XIS 1 detector. The 
first group of parameters are for data that include only statistical errors, while the second group of parameters are for data that also include 1.35% 
systematic errors. 


remain free. 

The results of these fits are presented in Fig. 5, and the pa- 
rameters are given in Table 2. Overall, the quality of the fits 
is very good, yielding 1756/1014= 1.74. Although 

the reduced is not 1 , the addition of only 1 .35% systematic 
uncertainties to each dataset brings the value down to that 
level, without significantly altering the best fit parameters, 
as also shown in Table 2. The addition of systematic un- 
certainties slightly increases the parameter error bars. 

Given the obvious systematic differences between the 
frontside (XIS 0 and 3) and backside (XIS 1) detectors, seen 
in the residuals presented in Fig. 5, we speculate that our fits 
are near “optimal’. Throughout the rest of this work, the 
figures will include only statistical errors; however, tables 
will include sets of parameters with error bars that were 
calculated for spectra to which 135% systematic errors 
have been applied. This systematic error amplitude was 
chosen on the basis that it yields a reduced value of « 1 
for the fits described above, and we retain this level in sub- 
sequent model fits for ease of comparison. 

The primary difference among the three spectra are the disk 
peak temperatures, which range between 1 .31-1 .48keV. The 
coronal contribution is relatively weak, given the low values 
of the coronal compactness compared to the disk compact- 
ness, 4/4- It is interesting to note that the two fainter ob- 
servations indicate less energetic coronae — values are 
lower — but higher optical depth ones. (Coronal compact- 
ness is its input enei^y divided by a characteristic radius, 
here measured relative to a compactness for the disk. If 
geometric changes are absent, then a lower fitted com- 
pactness is indicating a less enei^etic corona.) The result- 
ing equilibrium coronal temperatures for the three observa- 
tions are « 600 keV, « 2.2 keV, and ^ 1 .3 keV, respectively. 
When combined with their equilibrium optical depths, we find 
that the resulting Compton y-parameters range from 0. 1-0.2. 
That is, Comptonization is only inducing a modest shift in the 
average photon energy. However, given the fact that we do 
not have a good spectrum above 8 keV, there is not a strong 
enough lever arm on these Comptonization fits to say how 
free from systematic uncertainty these parameter values are. 

Likewise, detections of a (broad) Gaussian line peaking 
between 6.2-6.9keV are marginal. The upper limits of the 
line strengths found here (in terms of integrated photon flux) 
are comparable to the estimated values from our prior RXTE 
observations, after taking into account the fact that RXTE 


also integrates line emission from the galactic ridge and 
likely had remaining systematic Fe line residuals at the 1 % 
level (Nowak et al. 2008). That is, these line strengths are 
marginally consistent with the RXTE results. Furthermore, 
the line fitted here may be partly substituting for general spec- 
tral hardening at > 6 keV. It is clear, however, that any real line 
at the levels suggested by the RXTE observations cannot be 
narrow; any such line must be broad and have a peak residual 
relative to the continuum of < 2%. We illustrate this in Fig. 6, 
wherein we have refit the spectra excluding the Fe line, and 
combine all three observations in the plot — but not the fit — 
to show the Fe region residuals. 

The most salient feature of these fits is that a model with 
a single normalization constant is both successful, and yields 
a fairly small normalization of « 1.9 x 10^. The eqpair 
model normalization is proportional to the disk flux that 
would be observed in the absence of Comptonization, and is 
the same as the normalization of the diskpn model (Gier- 
lihski et al. 1999), namely 




eqp - 


M 


D 

1 kpc 




( 1 ) 


where M is the central object mass, D is the distance, i is the 
disk inclination, and /c is the color-correction factor of the 
disk (see also Coppi 1999) . For a mass of 3 Mq , a distance of 
10 kpc, and an inclination of i = 75°, we expect a normaliza- 
tion of fs 3 X 10“^, a factor of 15 times higher than observed. 
(The discrepancy only increases for a more face on geometry; 
a more highly inclined geometry is unlikely given the lack 
of eclipses.) The fitted normalization can be recovered if we 
increase the distance to 40 kpc, increase the color-correction 
factor for the seed input photons to 3.3, or we decrease the 
characteristic emission radius by a factor of four. (The emit- 
ting area dependence is responsible for the term in the 
normalization.) 

The diskpn model, which provides the seed photons 
for the eqpair model, has its peak emission radius at w 
\QGM/c^. If this radius were shrunk to 2.5GMjc^, the fitted 
normalization can be recovered with M = 3Mq, D = 10 kpc, 
and i = 75° . This is what we essentially now do by considering 
model fits with the kerrbb model. The kerrbb model (Li 
et al. 2005) uses a relativistic model of the disk structure, and 
includes Doppler beaming and gravitational redshift in calcu- 
lation of the emergent spectrum. The spin of the black hole 
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Fig. 7.— The joint fit of a simpl<g)kerrbb model to the three spectra 
assuming a black hole mass and distance of 16Mq and 22kpc. The spec- 
tra share a common interstellar absorption column, black hole mass and spin, 
disk inclination, and a fixed color correction factor of /c = 1 .7. They have 
independent disk accretion rates and simpl parameters. Black/grey spec- 
tra are the combined XIS 0 and 3 spectra, and blue/light blue are the XIS 1 
spectra. Lightly shaded data were not included in the fit. The top residuals 
panel shows the ratio residuals for obsID 405057010, the middle residuals 
panel shows the ratio residuals for obsID 405057020, and the bottom shows 
the ratio residuals for obsID 405057030. 

becomes a fit parameter. Increasing the spin from 0 (i.e., a 
Schwarzschild black hole) has two major effects: it increases 
the accretion efficiency of the disk, and decreases the disk in- 
ner radius. These two effects in combination tend to increase 
the disk temperature while decreasing its normalization (i.e., 
emitting area or, equivalently, inner radius). Both of these 
effects are suggested in the 4U 1957+1 1 spectra. 

Here we adopt the exact same set of parameters as we 
considered for our prior studies of the RXTE spectra of 
4U 1957+11 (Nowak et al. 2008). We assume a disk incli- 
nation of i = 75° , the theoretically preferred color-correction 
value of /c = 1.7 (Davis et al. 2005, 2006) for the disk spec- 
trum, and choose two sets of parameters: 3Mq and lOkpc, 
and 16M© and 22kpc. We discuss these latter two choice 
more fully in Section 4. We further modify these fits by al- 
lowing for a spectral hardening which we model with simpl 
(Steiner et al. 2009), which is a convolution model which 
mimics some aspects of optically thin Comptonization^. The 
two parameters of this convolution model are the slope of the 
powerlaw into which disk photons are transferred, and the 
fraction of input photons transferred into the hard tail. 

As for the Comptonization models, we fit all three observa- 
tions simultaneously, include a broad Gaussian line in exactly 
the same manner, and tie some of the model parameters to- 
gether. Specifically, we constrain each observation to be fit 
with the same neutral column, line energy and width (fixed 
to 0.3 keV), and black hole spin, but let the disk accretion 
rate, the simpl model parameters, and the line normaliza- 
tion be free. Results for these kerrbb model fits are pre- 


'' simpl, however, should not be confused with a Comptonization model. 
Aside from lacking any physical paramelere, e.g., a coronal optical depth, a 
coronal temperature or compactness, and any reference to coronal geometry, 
it also cannot mimic spectral shape.s associated with optically thick Comp- 
tonization. It should be regarded as purely a phenomenological model, albeit 
one that is more useful than a powerlaw since it does not extend its spectrum 
to low energies in the same manner as a steep powerlaw can. 


sented in Table 3 and Fig. 7. The resulting fits are good, with 
the 3Mq, lOkpc model yielding = 2063/1014 and the 
16M0, 22kpc model yielding y^^!v= 1830/1014. The latter 
model is very comparable to our best Comptonization fit. Ad- 
ditionally, this latter, more successful, model has line param- 
eters very similar to those from the Comptonization fits. We 
also find that only » 7 - 15% of the disk photons are trans- 
ferred into a hard tail. The fitted neutral column is slightly 
smaller than for the Comptonization fits, but only by < 10%. 

Just as the Comptonization models could fit all three obser- 
vations with a constant normalization (equivalent to a constant 
inner disk radius), so too can the relativistic disk models fit all 
three observations with a constant spin parameter (also equiv- 
alent to a constant inner disk radius). The main variable pa- 
rameter among the observations here becomes the disk accre- 
tion rate. For both sets of mass and distance assumptions, the 
fitted spin parameter is high: nearly 0.9 for the 3M©, lOkpc 
case, and nearly maximally spinning in the other. We discuss 
the scaling relations among the disk parameters in the next 
section. 


4. SCALING RELATIONS 

As shown above, spectral fits to the data are somewhat de- 
generate. We have found a successful set of relativistic disk 
model fits assuming a mass of 3 M© and a distance of 10 kpc, 
and also when assuming a mass of 16M© and a distance of 
22 kpc. The first set of parameters was chosen as being the 
lowest mass/closest distance black hole that is likely to be 
consistent with 4U 1957+1 1 persistently remaining in a soft 
state. The second set of parameters was chosen to be consis- 
tent with our prior RXTE studies (Nowak et al. 2008), which 
in turn were found via the scaling relations described below. 

As the spectra are dominated by the seed photons from a 
disk for both the disk and Comptonization models, the flux, 
F, scales as: 


F (X —T^f X 




D2 


ZkV 


D^KfcJ 


(2) 


where D is the distance to the source, R is the inner radius of 
the disk, M is the black hole mass, is the disk’s peak ef- 
fective temperature, T/oi is the disk’s peak color temperature, 
and /c is the color correction factor. The effective temper- 
ature to the fourth power scales as the fractional Eddington 
luminosity, L/Ledd. multiplied by M/R^ x M“'. Thus, if we 
make the assumption that the lowest flux RXTE observations 
described by Nowak et al. (2008) are indeed at a fixed fraction 
of Eddington luminosity, e.g., 3%, then we will have: 


Tool 


4 
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(3) 


The color temperature is fixed by observation, and L/Ledd « 
M/M (where M is the accretion rate) is fixed by assumption, 
leading to the first two of our scaling relations: 

Mxf, , Mxft . (4) 

Just as the color temperature is fixed by observations, so too 
is the flux. Given that the disk inner radius scales as mass, and 
using the fact that the color temperature and flux are fixed by 
observation, we have that M~ /EFf/ is fixed. Using the above 
scaling relations, we then obtain: 

Dxfl . 


(5) 
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TABLES 

constant*TBnew*simpl0{)cerrbb-fgaussian) Fits TO 4U 1957 -i-ll Spectra 


ObsID 

Ml 

a” 

M 

r 

/scat 


GxisO 

CxisS 

X+DoF 


(1(F‘ cm"-) 

(Jc/GM-) 




(10-*) 
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0.8965!«;““ 

0.180y«;“j 
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U.UUU_QQQQ 
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+ ^- 0.002 

1 O70+O-002 
i -u^U_o _002 
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j ,-K).] 
^•^- 0.0 
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1 A 70 +O.OO 3 
1 .U 003 

1 ART+0.003 
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0 1 1 5 + 0-002 
1 ->_ 0.001 


0015*^-^ 

^•V?IJ_OX)08 

0 O+O- 
^•^-0.0 

1 a^t+0.005 

l.UO/_Qyo_^ 

1 AA7+0.004 
1 ^-0.(K)5 
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n i«a+o.(X )2 

0.8957!)];“>li 

A 1 Rl+O.OOi 

, 1+2.9 
^•*- 0.0 

0-0005SS 

1 0+0.5 


1 qoi-H).003 

i-y-l_O.004 

1189.8/1014 

405057020 



U.l3i_0 002 

^•^-0.0 

0-032!«S 

0 

'^•'^-0.0 

1 077+^-^ 
-0.004 

1 AQ 1 -1-0.004 

l.Usi_o.(X)4 
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0.ii6tS5TT 

4.0!«“ 

0.009!«;»g 

0 0^'-^ 
^'^-0.0 

1 /Vi7+0-006 
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1 A^-i-0.006 

i -y^-o.ooo 
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A QQQ 1+0.0002 
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U.U /-_() (X)7 
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1 02 
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1 A-T 1+0-003 
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1830.0/1014 
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0 I^^)+O 007 
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1 
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1 .GOy_o,0()4 

1 n^R+O.004 
1 .U3o_0 005 


405057010 

p 

88 

0 QQ93+^ <^2 
^•^^■^- 0.0002 

0,282!«;0® 

1 4+0.2 

A A<!;»+ 0.010 

2.1T. 

1 A97+0.0&4 
1 -U— -0.003 
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0.219!»;“2 
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0.2^-^ 

1 Aar+0.(X)4 
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3 7+0.7 

0.128!«;“« 

A 9 + 0.6 
'^'^- 0.6 
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Note. — Errors are 90% confidence level for one interesting parameter. Ml is the neutral column; a* is the black hole dimensionless 
spin; M is the disk accretion rate; /c, the disk color-correction factor, has been fixed to 1.7; T and /scat are the simpl parameters 
(powerlaw slope, constrained to lie between 1.1-4, and fraction of scattered continuum); N-y is the Gaussian line normalization in units of 
photons cm”“ s“' ; and QisN are the normalization constants for the XIS 0 and 3 detectors, relative to the XIS 1 detector. The top set of 
parameters are for a black hole mass of SA/q and a distance of lOkpc, while the bottom set are for a mass and distance of 16Mq and 22kpc. 
Within each set of parameters, the first group are for data that include only statistical errors, while the second group are for data 
that also include 1.35% systematic errors. 


We have previously used these scaling relations to find 
spectral fits for a black hole mass of 16 M© and a distance of 

22kpc (Nowak et al. 2008). Specifically, we first explored 
fits to RXTE spectra where we fixed the mass and distance to 
3 M© and 10 kpc. We found that the best fits clustered around 
a spectral hardening factor of w 1 . 1 . We then fixed the spec- 
tral hardening factor to /c = 1.7 and the mass to 16M©. We 
allowed the distance to be a free parameter, and found that its 
best fit values clustered around 22 kpc, i.e., roughly consistent 
with the scaling relations above. 

Here we perform a similar exercise for the Suzaku spectra. 
As in Section 3 , we fit all three spectra simultaneously with 
the simpl^kerrbb model wherein we only allow the ac- 
cretion rate and simpl parameters to vary among the three 
spectra. We consider two sets of fits, one with mass of 3 Af© 
and distance 10 kpc, and another with 16 M© and distance 
22 kpc. For both we allow the black hole spin to be a free 
parameter, and we create grids of fits for fixed spectral hard- 
ening factors ranging over/c = 1 .1-3.7. The plot of fitted spin, 
a* , vs. spectral hardening factor, fc , is presented in Fig. 8. 

As for our prior RXTE observations, we find that fits us- 
ing 3M©, 10 kpc, and /c = 1.1 or 16M©, 22kpc,and /c = 1.7 
are nearly equally good and are of comparable to our best 
Comptonization fits . We also see that all fits using 1 6 M© , and 
22 kpc require large spin parameters. The fits using 3Mq and 
10 kpc only achieve low spin parameters for large color cor- 
rection factors, /c > 3. (These conclusions are unaltered if 
uniform 135% systematic errors are included in the fits.) 

5. DISCUSSION 

There are four salient features of these Suzaku observations 
of 4U 1957+1 1 that echo the results that we obtained with our 
prior Chandra, XMM-Newton, and RXTE studies. The peak 
disk temperatures are high, the disk normalizations (or equiv- 
alently, emitting areas or inner disk radii) are small, all ob- 
servations can be fit with the same inner disk radius, and the 
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Fig. 8.— Results of fitting a simpliSkerrbb to the three Suzaku ob- 
servations of 4U 1957+11. The top panel shows the fitted spin vs. spectral 
hardening factor. The bottom panel shows the reduced X" value of the fits. 
The blue (thicker) lines represent models with the mass and distance fixed to 
16 Mq and 22 kpc, respectively, and orange lines represent models with the 
mass and distance fixed to 3 Mq and 10 kpc, respectively. 

spectra are remarkably simple. The only significant devia- 
tions from a simple continuum, aside from spectral hardening 
and a possible weak, broad Fe line at high energy, are « 0.7- 
0.9 keV features associated with absorption by the ISM, pre- 
dominantly due to the Fe L and Ne edges. 

We show a closeup of this region in Fig. 9. The spec- 
tral feature seen at sa 0.55 keV, not described by our fitted 
models, can in fact be seen in the flux corrected spectra of 
any monochromatic, narrow line spectrum simulated using 
the Suzaku response functions. Small calibration errors in 
the amplitude of this feature potentially could lead to a sys- 
tematic underestimate or overestimate of the predicted counts 
near 0.55 keV. We thus do not believe that this feature is in- 
dicative of physical complexity in the underlying 4U 1957+11 
spectra. Likewise, it is unclear whether the w 5% residuals 
seen near the 0.7 keV Fe L edge represent complexity in the 
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FlO. 9.— The same fit as seen in Fig. 5, but showing the XIS 0, 1, and 3 
spectra combined, and the model applied down to 0.5 keV. 


If one were to reduce this radius to 2.5 GM/(p-, i.e., change 
the emission profile to more closely match that of a rapidly 
spinning hole, then the Comptonization model normalization 
would be consistent with a 3 Mq black hole at only 10 kpc. 

The Comptonization models offer two possibilities® for fit- 
ting the 4U 1957 -fI 1 spectrum with a low spin. The first, per- 
haps less likely, possibility is to increase the color-correction 
factor via a low temperature, optically thick corona. We have 
already seen that the lowest flux observations are fit with such 
coronae, albeit with a Compton y-parameter of « 0. 1 . Instead, 
we are envisioning a scenario wherein y « /c « 3 - the color 
correction value required to fit a 3 Mq , 1 0 kpc, Schwarzschild 
black hole. This would require an optical depth of 


r 


9.8 


TkeVV''^ 



(6) 


physical spectrum, as opposed to uncertainty in the response. 
Aside from the Fe L edges, Chandra and XMM-Newton spec- 
tra did not show any spectral features near the 0.7 keV band- 
pass (Nowak et al . 2008). 

The only clear spectral features are the Ne edge/Ne IX ab- 
sorption near 0.9 keV. Similarly simple disk spectra were im- 
plied by Suzaku observations of LMC X-3 in its soft, disk 
dominated state (Kubota et al. 2010). As pointed out by Kub- 
ota et al. (2010), such spectral simplicity is in fact surpris- 
ing if modem disk atmosphere models apply to these sources. 
These disk atmosphere models predict spectral structure with 
±5% amplitude near metal edges (Davis & Hubeny 2006) that 
is seen neither here nor in Suzaku observations of LMC X- 
3. That the simplest, phenomenological disk models lack- 
ing such features describe the spectra extremely well (i.e., to 
within < 5% 0 .7-0 .8 keV, and to < 2% above 0.8 keV) is quite 
remarkable, and may indicate that additional physics needs 
to be incorporated into the disk atmosphere models (Kubota 
etal. 2010). 

We also have confirmed the results of our prior studies in 
that the best fits assuming a small mass and distance (3 Mq, 
10 kpc) require a black hole spin a* > 0.9. The required spin 
increases if we instead assume 16Mq and 22 kpc. Altering 
the assumed disk inclination does not change this situation. 
Higher disk inclinations are ruled out by the lack of system 
eclipses. Lower disk inclinations allow for more observed 
gravitational redshifting, and less relativistic beaming, of 
the spectrum, and thus require larger spectral correction fac- 
tors. (Le., inclination effects more than just the projected 
area of the disk.) As an example, we took the 16Mq and 
22 kpc kerrbb model, froze the inclination to 45° and let 
both the black hole spin and hardness correction factor be- 
come free parameters. A solution with comparable to 
the best fits discussed above was found, however, it required 
maximally negative spin of a* ^ -1 and an extremely large 
color correction factor of « 8. Thus, high spin/high incli- 
nation solutions are the most “straightforward” solutions that 
we have been able to obtain, and generally have had the 
lowest color-correction factors of the models that we have 
explored (see also Nowak et al. 2008). 

These results are in fact confirmed by the Comptoniza- 
tion model fits: the high disk temperatures (1.31-1 .48keV) 
and especially the low normalization — 15 times lower than 
expected for a SM© Schwarzschild black hole at 10 kpc — 
could be indicative of a rapidly spinning black hole. As dis- 
cussed above, the disk temperature in the diskpn seed pho- 
tons used in the Comptonization model peaks at w 10 GM/c^. 


where kT^ is the coronal temperature. This is a coronal regime 
typically not explored in spectral models. The concept would 
be that the low flux, equilibrium solution of the disk was such 
a cool, optically thick atmosphere which at higher flux (i.e., 
the 15% of the time requiring a hard tail) expanded into a hot, 
optically thin corona. 

Aside from lacking a theoretical description of such a sce- 
nario, the fitted constant emitting area at low flux implies that 
the 4U 1957 -fI 1 system must always return to the same equi- 
librium coronal solution at low flux. It is not clear why this 
should be the case. The second possibility offered by the 
Comptonization solutions is simply to increase the distance 
to 4U 1957+1 1 , while keeping the mass small. The scaling 
relations discussed in Section 4 only apply under the assump- 
tion that the lowest flux observations are a fixed fraction of 
the Eddington luminosity (e.g., 2%LEdd)- If we relax this 
assumption, then eq. 1 suggests that if the 4U 1957+1 1 dis- 
tance were 40 kpc, then it could be a 3M© Schwarzschild 
black hole. This would imply, however, that the faintest 
4U 1957+1 1 spectra are > 30%LEdd> while the brightest ob- 
served RXTE spectra are > LEdd- This would be rather un- 
usual, especially given the lack of strong variability found in 
the RXTE spectra: the root mean square variability is typi- 
cally < 3%, and at its strongest was < 13% in the 2-60 keV 
bandpass over « 10'^-2 x 10^ Hz (Nowak & Wilms 1999; 
Nowak et al. 2008). 

It therefore seems likely that 4U 1957+1 1 is a system that 
has a combination of both high spin and large distance, which 
puts this X-ray binary in the Galaxy’s halo. It is interesting 
to speculate whether these two characteristics are related to 
one another. First, we note the obvious fact that in order for 
4U 1957+11 to be a X-ray binary it must still be accreting 
and therefore be a comparatively young system. Given the 
low star formation rate in the Galaxy’s halo, it is rather un- 
likely, therefore, that 4U 1957+11 was formed in the halo. 
Instead 4U 1957+1 1 likely originates in the Galactic disk and 
then migrated into the halo. Several possible formation sce- 
narios have been discussed in the literature, most recently 
associated with the discovery of young hyperv'elocity stars 
in the Galaxy’s halo. For example, Perets (2009, and refer- 
ences therein) discuss a scenario where hypervelocity binary 
systems are created through the interaction of a hierarchical 

® Note that we are ignoring issues of the unknown inclination. Choosing 
a more face-on inclination only exacerbates the normalization problem (i.e., 
we would have expected an even larger normalization), and would introduce 
more gravitational redshifting to the observed temperature, increasing the re- 
quirement for rapid spin in the spectral fits. 
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triple star system with the supennassive black hole in the cen- 
ter of the Galaxy. Since hierarchical triples are very com- 
mon among multiple star systems, such a scenario is likely. 
With a typical ejection velocity of ISOOkms"', it would take 
this system a few 10* years to reach its current position, well 
within the lifetime of typical evolution scenarios for X-ray 
binary systems (van den Heuvel 1976). (At such a veloc- 
ity, and assuming a 22 kpc distance, the proper motion of 
4U 1957+11 likely would have amounted to ^ 0.3" since 
the discovery of the optical counterpart by Margon et al. 
1978.) Such a scenario could also explain a large spin for 
the system, if during the ejection angular momentum is trans- 
ferred onto the black hole. 

A high black hole spin for4U 1957+1 1 is a consistent inter- 
pretation of the results presented here; however, such a result 
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